Measurement of the proton-air cross-section at y/s = 57TeV 
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'*Los Alamos National Laboratory, Los Alamos, NM, USA 
(Dated: August 10, 2012) 

We report a measurement of the proton-air cross-section for particle production at the center-of- 
mass energy per nucleon of 57TeV. This is derived from the distribution of the depths of shower 
maxima observed with the Pierre Auger Observatory: systematic uncertainties are studied in de- 
tail. Analysing the tail of the distribution of the shower maxima, a proton-air cross-section of 
[505 ± 22(stat) ^(sys)] mb is found. 

PACS numbers: 



INTRODUCTION 

We present the first analysis of the proton-air cross- 
section based on measurements made at the Pierre Auger 
Observatory [l[. For this purpose we analyse the shape 
of the distribution of the largest values of the depth of 
shower maximum, X max , the position at which an air 
shower deposits the maximum energy per unit of mass of 
atmosphere traversed. The tail of the X max -distribution 
is sensitive to the proton-air cross-section, a fact first 
exploited in the pioneering work of the Fly's Eye Collab- 
oration 0. To obtain accurate measurements of X max , 
timing data from the fluorescence telescopes is combined 
with that from the surface detector array for aprecise 
hybrid reconstruction of the geometry of events |3[. 

We place particular emphasis on studying systematic 
uncertainties in the cross-section analysis. The unknown 
mass composition of cosmic-rays Q is identified to be 
the major source of systematic uncertainty and accord- 
ingly the analysis has been optimised to minimise the 
impact of particles other than protons in the primary 
beam. This begins with restricting the analysis to the 
energy interval 10 18 to 10 185 eV, where the shape of the 
X max -distribution is compatible with there being a sub- 
stantial fraction of protons; also there are a large number 
of events recorded in this energy range. The correspond- 
ing average center-of-mass energy of a proton interacting 
with a nucleon is 57TeV, significantly above the reach of 
the Large Hadron Collider. 

ANALYSIS APPROACH 

The proton-air cross-section is derived in a two-step 
process. Firstly, we measure an air-shower observable 
with high sensitivity to the cross-section. Secondly, we 
convert this measurement to a value of the proton-air 
cross-section for particle production (c.f. Q). This is 
the cross-section that accounts for all interactions which 
produce particles and thus contribute to the air-shower 
development; it implicitly also includes diffractive inter- 
actions. As the primary observable we define A,, via the 
exponential shape of the tail of the AT max -distribution, 
dN/dX mliK cx exp(— X ma , x / A v ), where r\ denotes the frac- 
tion of most deeply penetrating air showers used. Con- 
sidering only these events enhances the contribution of 



protons in the sample, since the depth at which proton- 
induced showers maximise is deeper in the atmosphere 
than for showers from heavier nuclei. Thus, r\ is a key 
parameter: a small value enhances the proton fraction, 
but reduces the number of events available for the analy- 
sis. We have chosen r\ = 0.2 so that, for helium-fractions 
up to 25 %, biases introduced by the possible presence of 
helium and heavier nuclei do not exceed the level of the 
statistical uncertainty. This was chosen after a Monte 
Carlo study that probed the sensitivity of the analysis to 
the mass composition depending on the choice of different 
values of r/. 

THE MEASUREMENT OF A,, 

We use events collected between 1 Dec 2004 and 20 
Sept 2010. The atmospheric and event-quality cuts ap- 
plied are identical to those used for the analysis of (X max ) 
and RMS(A max ) Q yielding 11628 high-quality events. 
The AT max -distribution of these data is affected by the 
known geometrical acceptance of the fluorescence tele- 
scopes as well as by limitations related to atmospheric 
light transmission. We use the strategy developed for 
the measurement of (AT max ) and RMS(A" max ) to extract 
a sample that has an unbiased X max -distribution: a fidu- 
cial volume selection, which requires event geometries 
that allow, for each individual shower, the complete ob- 
servation of a defined slant depth range. 

Firstly, we derive the range of values of X max that 
corresponds to the fraction T) — 0.2 of the most deeply 
penetrating showers. For this we need an unbiased distri- 
bution of A" max over the entire depth range of observed 
values of A" max . To achieve this we perform a fiducial 
event selection of the slant depth range containing 99.8 % 
of the observed X max -distribution, which corresponds to 
the range from 550 to 1004 g/cm 2 . This reduces the 
data sample to 1635 events providing an unbiased A" max - 
distribution that is used to find the range of values of 
I max corresponding to T) = 0.2, identified to extend from 
768 to 1004 g/cm 2 . 

Secondly, we select those events from the original 
11628 that have geometries allowing the complete ob- 
servation of values of AT max from 768 to 1004 g/cm 2 , the 
tail of the unbiased distribution. This fiducial cut max- 
imises the statistics of an unbiased AT max -distribution in 
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FIG. 1: Unbinned likelihood fit to obtain A,, (thick line). 
The Xmax-distribution is unbiased by the fiducial geometry 
selection applied in the range of the fit. 



the range of interest. In total 3082 events pass the fidu- 
cial volume cuts, of which 783 events have their A" max 
in the selected range and thus contribute directly to the 
measurement of A,,. In Fig. [I] we show the 3082 selected 
events and the result of an unbinned maximum likeli- 
hood fit of an exponential function over the range 768 to 
1004 g/cm 2 . Values of A^ have been re-calculated for sub- 
samples of the full dataset selected according to zenith- 
angle, shower-to-telescope distance and energy: the dif- 
ferent values obtained for A^ are consistent with statisti- 
cal fluctuations. The re-analyses of the data for changes 
of fiducial event selection, modified values of r\ and for 
different ranges of atmospheric depths yield changes of 
A,, that are distributed around zero with a root-mean- 
square of 1.6 g/cm 2 . We use this root-mean-square as an 
estimate of the systematic uncertainties associated to the 
measurement. This yields 



A r) = 



[55.8 ± 2.3(stat) ± 1.6(sys)] g/cm 2 



(1) 



with the average energy of these events being 
10 is.24 ±o.oo5(stat) cV ^ The differential energy distribution 

for these events follows a power-law with index —1.9. The 
average energy corresponds to a center-of-mass energy of 
-y/s = [57 ± 0.3(stat)] TeV in proton-proton collisions. 



DETERMINATION OF THE CROSS-SECTION 

The determination of the proton-air cross-section for 
particle production requires the use of air-shower sim- 
ulations, which inherently introduces some dependence 
on model assumptions. We emulate the measurement of 
A^ with Monte Carlo simulations to derive predictions of 
the slope, A^ c . It is known from previous work that the 
values of Aj^ c are directly linked to the hadronic cross- 
sections used in the simulations 0] . Accordingly we can 



explore the effect of changing cross-sections empirically 
by multiplying all hadronic cross-sections input to the 
simulations by an energy-dependent factor Q 



f(E, f 19 ) = 1 + (/i 9 - 1) 



lg(£/10 15 eV) 
lg(10 19 cV/10 15 eV)' 



(2) 



where E denotes the shower energy and /ig is the factor 
by which the cross-section is rescaled at 10 19 eV. This 
factor is unity below 10 15 eV reflecting the fact that mea- 
surements of the cross-section at the Tevatron were used 
to tune the interaction models. This technique of modi- 
fying the original predictions of the cross-sections during 
the simulation process assures a smooth transition from 
accelerator data up to the energies of our analysis. 

For each hadronic interaction model, the value of /ig is 
obtained that reproduces the measured value of A^ . The 
modified cross-section is then deduced by multiplying the 
original cross-section used in the model by the factor 
f(E, Jig) of Eq. @ using E = 10 18 24 eV. For the conver- 
sion of A, ; into cross-section, we have used the four high- 
energy hadronic interaction models commonly adopted 
for air shower simulations: QGSJetOl H, QGSJetII.3 [§], 
SIBYLL 2.1 and EPOS1.99 [ll[. While in general 
no model gives a completely accurate representation of 
cosmic-ray data in all respects, these have been found to 
give reasonably good descriptions of many of the main 
features. It has been shown [l2| that the differences be- 
tween the models used in the analysis are typically bigger 
than the variations obtained within one model by param- 
eter variation. Therefore we use the model differences for 
estimating the systematic model dependence. 

The proton-air cross-sections for particle production 
derived for QGSJetOl, QGSJctll, SIBYLL and EPOS 
are 523.7, 502.9, 496.7 and 497.7 mb respectively, with 
the statistical uncertainty for each of these values being 
22 mb. The difference of these cross-sections from the 
original model predictions arc < 5%, with the exception 
of the result obtained with the SIBYLL model, which is 
12 % smaller than the original SIBYLL prediction. Wc 
use the maximum deviations derived from using the four 
models, relative to the average result of 505 mb, to es- 
timate a systematic uncertainty of (—8, +19) mb related 
to the difficulties of modelling high energy interactions. 
This procedure relies on the coverage of the underlying 
theoretical uncertainties by the available models. For 
example diffraction, fragmentation, inelastic intermedi- 
ate states, nuclear effects, QCD saturation, etc. are all 
described at different levels using different phenomeno- 
logical, but self-consistent, approaches in these models. 
It is thus possible that the true range of the uncertainties 
for air-shower analyses is larger, but this cannot be esti- 
mated with these models. Furthermore, certain features 
of hadronic particle production, such as the multiplicity, 
elasticity and pion-chargc ratio, have an especially im- 
portant impact on air shower development (l3l . 14 1: of 



6 



related to the conversion of A,, to cr^^.. 

The presence of photons in the primary beam would 
bias the measurement. The average X max of showers 
produced by photons at the energies of interest is about 
50 g/cm 2 deeper in the atmosphere than that of protons. 
However, observational limits on the fraction of photons 
are < 0.5% [l6, 17 1. With simulations we find that the 
possible under-estimation of the cross-section if photons 
were present in the data sample at this level is less than 
10 mb. 

With the present limitations of observations, we can- 
not distinguish air showers produced by helium nuclei 
from those created by protons. From simulations we find 



TABLE I: Summary of the systematic uncertainties. 



Description 


Impact on o^_°^ T 


A,, systematics 


±15 mb 


Hadronic interaction models 




Energy scale 


±7mb 


Conversion of A,, to cr^°^ 


±7mb 


Photons, <0.5% 


< +10 mb 


Helium, 10 % 


-12 mb 


Helium, 25 % 


-30 mb 


Helium, 50 % 


-80 mb 


Total (25 % helium) 


-36 mb, +28 mb 



500 



400 



300 



200 V 
10" 



these we found that only the elasticity can have a rele- 
vant impact on A^. The previously identified systematic 
uncertainty of (—8, +19) mb induced by the modelling of m 700 
hadronic interactions, corresponds to the impact of mod- £ 
ifying the elasticity within ±(10 — 25) % in the models. c- 600 

The selection of events with large values of X max also § 

enhances the fraction of primary cosmic-ray interactions i, 

with smaller multiplicities and larger elasticities, which | 

is for example characteristic for diffractive interactions. I 

The value of A„ is thus more sensitive to the cross-section » 

o 

of those interactions. The identified mo del- dependence o 
for the determination of <Jp™t T i s a l so caused by the com- 
pensation of this effect. 

Also the choice of a logarithmic energy dependence 
for the rescaling-factor in Eq. @ may affect the result- 
ing cross-sections. However, since the required rescaling- 
factors are small, this can only be a marginal effect. 

The systematic uncertainty of 22 % [l5[ in the absolute 
value of the energy scale leads to systematic uncertainties 
of 7 mb in the cross-section and 6 TeV in the center-of- 
mass energy. Furthermore, the procedure to obtain u^™? 
from the measured A,, depends on additional parameters. 
By varying the energy distribution, energy and AT max res- 
olution in the simulations, we find that related system- 
atic changes of the value of Cp™ ir are distributed with a 
root-mean-square of 7mb around zero. We use the root- 
mean-square as estimate of the systematic uncertainties 



Equivalent cm. energy is pp [TeV] 
Iff 1 1 10 



Nam et al. 1975 [30] 
Siohan el al. 1978 [31] 
Baltrusaitis el al. 1 984 [2] 
Mielkeetal. 1994 [32] 
Knurenko et al. 1999 [19] 
Honda el al. 1999 [20] 
Belov et al. 2007 [18] 
Aglietta et al. 2009 [33] 
Aielli et al. 2009 [34] 
■ This work 




QGSJetOlc 

QGSJetll.3 

Sibyll 2.1 

Epos 1.99 



10 1 ' 



10 1 " 



10 1 ' 



I I I 

10 15 10 16 10 1 "' 
Energy [eV] 



10" 



10" 



10^ 



FIG. 2: Resulting cr 



prod 



compared to other measure- 



ments (see [18H20I . I30H34I ] for references) and model predic- 
tions. The inner error bars are statistical, while the outer 
include systematic uncertainties for a helium fraction of 25 % 
and 10 mb for the systematic uncertainty attributed to the 
fraction of photons. 



that fJp™ ir is over-estimated depending on the percent- 
ages of helium in the data sample. Lack of knowledge of 
the helium fraction is the dominant source of systematic 
uncertainty. 

We also find that the nuclei of the CNO-group intro- 
duce no bias for fractions up to ~ 50 %, and accordingly 
we assign no uncertainty in the cross-section due to these 
or heavier nuclei. 

In Table H] we list the sources of systematic uncertain- 
ties. As the helium fraction is not known we show the 
impact of 10, 25 and 50 % of helium respectively. In what 
follows we include a systematic uncertainty related to a 
helium fraction of 25%. In the extreme case, were the 
cosmic-ray composition to be 100 % helium, the analysis 
would over-estimate the proton-air cross-section by 300 
to 500 mb. Given the constraints from accelerator data 
at lower energies and typical model assumptions, this ex- 
treme scenario is not realistic. 

We summarise our results by averaging the four values 
of the cross-section obtained with the hadronic interac- 
tion models to give 

- prod - [505 ±22(stat) t^(sys)] mb 



p-air 



at a center-of-mass energy of [ 57 ± 0.3(stat) ± 
6(sys) ] TeV. In Fig. [2] we compare this result with model 
predictions and other measurements. The measurements 
at the highest energies are: HiRes [HI and Fly's eye 
that are both based on X max , Yakutsk Array [l9| using 
Cherenkov observations and Akeno [2(| measuring elec- 
tron and muon numbers at ground level. All these analy- 
ses assume a pure proton composition. In the context of 
a possible mixed-mass cosmic-ray composition, this can 
lead to large systematic effects. Also all these analyses 
are based on a single interaction model for describing air 
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showers: Only HiRes uses a second model for systematic 
checks. 

It is one of the prime aims of our analysis to have the 
smallest possible sensitivity to a non-proton component, 
and to perform a detailed systematic analysis on the un- 
certainties related to the mass composition. We also use, 
for the first time, all hadronic interaction models cur- 
rently available for the estimation of model-related sys- 
tematic effects. Futhermorc, by using Eq. ^ty we derive 
a cross-section corresponding to a smooth interpolation 
from the Tevatron measurement to our analysis, with no 
inconsistencies as in earlier approaches. 



COMPARISON WITH ACCELERATOR DATA 

For the purpose of making comparisons with acceler- 
ator data we calculate the inelastic and total proton- 
proton cross-sections using the Glauber model. We 
use standard Glauber formalism [2l|, extended by a 
two-channel implementation of inelastic intermediate 
states @ to account for diffraction dissociation [22[ ■ The 
first channel corresponds to p — > p scattering and has 
an amplitude of T pp , while the amplitude for the other 
channel is T pp » = A T pp and corresponds to the excitation 
of a short lived intermediate state. The parameter A is 
related to the ratio of singlc-diffractive cross-section and 
clastic cross-section. We use a value of A = 0.5±0.15 that 
is determined from measurements of the single-diffractivc 
cross-section, as well as from proton-carbon cross-section 
data at lower energies. 

This Glauber calculation is model-dependent since nei- 
ther the parameters nor the physical processes involved 
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FIG. 3: Correlation of elastic slope parameter, B e \, and the 
inelastic proton-proton cross-section in the Glauber frame- 
work. The solid line indicates the parameter combinations 
yielding the observed proton-air production cross-section, and 
the dotted lines are the statistical uncertainties. The hatched 
area corresponds to the predictions by SIBYLL, QGSJet, 
QGSJetll and EPOS. See also Ref. [|. 



are known accurately at cosmic-ray energies. In particu- 
lar, this applies to the clastic slope parameter, B c \. (de- 
fined by da e i/dt oc exp(—\t\B e \) for very small t), the cor- 
relation of B c \ to the cross-section, and the cross-section 
for diffractive dissociation. For the example of <r pp cl , the 
correlation of B e \ with the cross-section is shown in Fig. [3] 
for A = 0.5. We have used the same four hadronic inter- 
action models to determine the uncertainty band of the 

B e l_ 

as 



r inel 

PP 



23] fall within this band. We find that in the Glauber 



framework the inelastic cross-section is less dependent 
on model assumptions than the total cross-section. The 
result for the inelastic proton-proton cross-section is 

a™ 1 = [92 ±7(stat) ±n(sys) ± 7(Glauber)] mb, 

and the total proton-proton cross-section is 

a p f = [133 ± 13(stat) ±H(sys) ± 16(Glaubcr)] mb. 

The systematic uncertainties for the inelastic and total 
cross-sections include contributions from the elastic slope 
parameter, from A, from the description of the nuclear 
density profile, and from cross-checking these effects us- 
ing QGSJetll For the inelastic case, these three 
independent contributions are 1, 3, 5, and 4mb respec- 
tively. For the total cross-section, they are 13, 6, 5, and 
4mb. We emphasize that the total theoretical uncer- 
tainty of converting the proton-air to a proton-proton 
cross-section may be larger than estimated here within 
the Glauber model. There are other extensions of the 
Glauber model to account for inelastic screening d, [25[ 
or nucleon- nucleon correlations (2(|, and alternative ap- 
proaches that include, for example, parton saturation or 
other effects [H IS H3 . 

In Fig. |4] we compare our inelastic cross-section result 
to accelerator data and to the cross-sections used in the 
hadronic interaction models. 



SUMMARY 

We have presented the first measurement of the cross- 
section for the production of particles in proton-air col- 
lisions from data collected at the Pierre Auger Observa- 
tory. We have studied in detail the effects of assumptions 
on the primary cosmic-ray mass composition, hadronic 
interaction models, simulation settings and the fiducial 
volume limits of the telescopes on the final result. By 
analysing only the most deeply penetrating events we se- 
lected a data sample enriched in protons. The results are 
presented assuming a maximum contamination of 25 % 
of helium in the light cosmic-ray mass component. The 
lack of knowledge of the helium component is the largest 
source of systematic uncertainty. However, for helium 
fractions up to 25 % the induced bias remains below 6 %. 

To derive a value of 0^'°^ from the measured A v we 
assume a smooth extrapolation of hadronic cross-sections 
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FIG. 4: Comparison of derived a™ 1 to model predictions and 
accelerator data [291 ]. Here we also show the cross-sections of 
two typical high-energy models, Pythia6 [3^] and Phojet[3(|. 
The inner error bars are statistical, while the outer include 
systematic uncertainties. 



from accelerator measurements to the energy of the anal- 
ysis. This is achieved by modifying the model-predictions 
of hadronic cross-sections above energies of 10 15 eV dur- 
ing the air-shower simulation process in a self-consistent 
approach. 

We convert the proton-air production cross-section 
into the total, and the inelastic, proton-proton cross- 
section using a Glauber calculation that includes inter- 
mediate inelastic screening corrections. In this calcula- 
tion we use the correlation between the elastic slope pa- 
rameter and the proton-proton cross-sections taken from 
the interaction models as a constraint. We find that the 
inelastic proton-proton cross-section depends less on the 
clastic slope parameter than docs the total proton-proton 
cross-section, and thus the systematic uncertainty of the 
Glauber calculation for the inelastic result is smaller. 
The data agree with an extrapolation from LHC [29| en- 
ergies to 57TeV for a limited set of models. 
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